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Abstract. We study numerically the monopole creation operator proposed 
recently by Frohlich and Marchetti. The operator is defined with the help 
of a three dimensional model which generates random Mandelstam strings. 
These strings imitate the Coulombic magnetic field around the monopole. 
We show that if the Mandelstam strings are condensed the creation oper- 
ator discriminates between the phases with condensed and non-condensed 
monopoles in the Abelian Higgs model with the compact gauge field. 
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1. Introduction 

The order parameters are very important for the investigation of the phase 
transitions. In case of the temperature phase transition in lattice QCD with 
dynamical quarks there is no good definition of the order parameter up to 
now. The traditional order parameters like Polyakov line and string tension 
which work well in the quenched case are not valid in case of the full QCD. 
Even at zero temperature the string between quark and anti-quark can be 
broken by sea quarks. 

We discuss below a quantity which may serve as the order parameter for 
full QCD if the monopole (or, "dual superconductor") confinement mecha- 
nism [1] is valid. In this picture the monopoles - defined with the help of an 
Abelian projection [2] - are supposed to be condensed in the confinement 
phase. The monopole condensate causes a dual analogue of the Abrikosov 
vortex to be formed between quarks and anti-quarks. As a result the quarks 
and anti-quarks are confined into the colorless states. In the deconfinement 
phase the monopoles are not condensed and quarks are not confined. Thus 
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the natural confinement-deconfinement order parameter is the value of the 
monopole condensate, which should be nonzero in the confinement phase 
and zero in the deconfinement phase. 

There are two (formal) difficulties in the definition of the monopole 
condensate. At first, the monopole condensate is defined as an expectation 
value of the monopole field. However, the monopoles are the topological 
defects in the compact Abelian gauge field and, as a result, the immediate 
output of the lattice simulations is an information about the monopole 
trajectories. Then one can apply a known procedure which allows us to 
rewrite the path integral over monopole trajectories as an integral over the 
monopole fields. The expectation value of the latter defines the monopole 
condensate. 

The next difficulty is the following. The expectation value of the (e.g., 
scalar) field <fi should always be zero regardless whether this field is con- 
densed or not. The reason is very simple: the path integral includes the 
integration over all possible gauges while the charged field is gauge de- 
pendent. These two problems were solved in paper [3], where the gauge 
invariant monopole creation operator for compact QED (cQED) was ex- 
plicitly constructed. The numerical calculations in lattice cQED [4] and in 
the Maximal Abelian projection of lattice ST/ (2) gluodynamics, Ref. [5], 
show that this operator provides the order parameter for confinement - 
deconfinement phase transition. Later the other forms of the monopole cre- 
ation operators were constructed and investigated numerically [6] in lattice 
gluodynamics. 

But recently it appears the claim [7] that the monopole creation opera- 
tor suggested in [3] depends on the shape of the Dirac string if dynamical 
electrically charged fields are present in the considered theory. Just this sit- 
uation appears in the Abelian projection of QCD: the off-diagonal gluons 
become electrically charged dynamical fields while the diagonal gluons be- 
come compact Abelian gauge fields which contain monopoles. The authors 
of [7] suggested a "new" monopole creation operator, which does not de- 
pend on the shape of the Dirac string even in the presence of the dynamical 
electric charges. 

Below we study this new creation operator in the compact Abelian Higgs 
model, having in mind the future application of this operator for the full 
QCD. 

In the next Section we give the explicit construction of the "new" and 
the "old" monopole creation operator. In Section 3 we present the results 
of the numerical calculations in the compact Abelian Higgs model. 
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2. "Old" and "new" monopole creation operators. 

The gauge invariant creation operator <E> was suggested by Dirac [8]: 

$ = <£(x)expji J E k (x - y)A k (y) d 3 y^ , (1) 

here 4>(x) and A k (x) are the electrically charged field and gauge potential, 
which transform under the gauge transformations as 

<f>(x) - <t>{x)e ia ^ , A k (x) ^ A k (x) + d k a(x) . (2) 

The Coulomb field, E k (x), satisfies the equation: 

d k E k = 5^(x). (3) 

It is easy to see that the operator $, eq. (1), is invariant under the gauge 
transformations (2). 

Now we describe the Frohlich-Marchetti construction [3] of the monopole 
creation operator in cQED. At first step the partition function of cQED is 
transformed to a dual representation. For the general form of the cQED 
action it can be shown [5] that the dual theory is an Abelian Higgs model 
(AHM) in the limit when the Higgs boson mass and the gauge boson mass 
are infinite. In this theory the Higgs field, <fi x , corresponds to the monopole 
in the original cQED. The gauge field *B is dual to the original gauge 
field 9. Thus the gauge invariant creation operator (1) for the AHM model, 
corresponds to the monopole creation operator in the original cQED. The 
explicit expression for this operator on the lattice is (cf. eq.(2)): 

= <t> x e i( -* B * H *\ (4) 

where *H X is the Coulomb field of the monopole, 5*H X = *5 X , and *5 X is the 
discrete 5-function defined on the dual lattice. Here and below we will use 
the differential form notations on the lattice: (a, b) = J2c a ct>c is the scalar 
product of the forms a and b defined on the c-sells; (a, a) = ||a|| 2 is the 
norm of the form a; d is the forward derivative (an analog of the gradient); 
5 is the backward derivative (an analog of the divergence) and ^-operation 
transfers a form to the dual lattice. For a description of the language of the 
differential forms on the lattice see, e.g., review [9]. 

Performing the inverse duality transformation for the expectation value 
of the creation operator (4) we get the expectation value, ($ mon ), of this 
operator in cQED. The explicit expression in the lattice notations is: 

($ mon ) = \ fvO exp{-S(d0 + W)} , 

Z = fve exp{-S(d0)}, (5) 

J — 7T 
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here d9 is the plaquette angle, the lattice action is a periodic function: 
S(d9 + 27m) = S(d9), n € Z\ W = 2k5A~ 1 *(H x - u x ) and *u x is the 
Dirac string which starts at the monopole: 5*lo x = *8 X . The Dirac string 
*lo x is defined on the dual lattice. The numerical investigation of this cre- 
ation operator in cQED shows [4] that it can be used as the confinement- 
deconfinement order parameter. 

The operator (4) is well defined for the theories without dynamical 
matter fields. However, if an electrically charged matter is added, then the 
creation operator (4) depends on the position of the Dirac string. To see 
this fact let us consider the compact Abelian Higgs model with the Villain 
form of the action: 

Zahm= fvO fvp E e -' 3 ll de + 2 ™ll 2 ^ll d ^ +M H 2 . (6) 

Here 9 is the compact Abelian gauge field and <p is the phase of the dy- 
namical Higgs field. The integer q is the charge of the Higgs field. For the 
sake of simplicity we consider the London limit (the Higgs mass is infinitely 
large while the Higgs condensate is constant). 

Let us perform the Berezinsky-Kosterlitz-Thouless (BKT) transforma- 
tion [10] with respect to the compact gauge field 9: 

d9 + 2vrn = dA + 2irSA~ 1 j , with A = 9 + 2vr5A" 1 m[j] + 2nk . (7) 

Here A is the non-compact gauge field, *m[j] is a surface on the dual 
lattice spanned on the monopole current *j (5*m\j] = *j), A is the lattice 
Laplacian and k is the integer-valued vector form 1 . We substitute eqs.(7) 
in eq.(6) and make the shift of the integer variable, I — > I + qk. 

Next we perform the BKT transformation with respect to the compact 
scalar field ip: 

dtp + 2vr/ = d-d + 2tt(5A" 1 ct , with ■& = ip + 27t5A~ 1 s[ct] + 2tt P . (8) 

Here t? is the non-compact scalar field, *s[a] is a 3D hyper-surface on the 
dual lattice spanned on the closed surface V (<5*s[<r] = V) and p is the 
integer-valued scalar form. 

Substituting eqs. (7,8) into the partition function (6) and integrating 
the fields A and p we get the representation of the compact AHM in terms 
of the monopoles and strings ("the BKT-representation" ) : 

Z AHM <xZ BKT = ]T ]T exp{ -4tt 2 /3 (j,(A + m 2 ) _1 i) (9) 
S*j=0 g* aj=q *j 

: A detailed description of the duality and BKT transformations in terms of the dif- 
ferential forms on the lattice can be found, e.g., in Ref. [9]. 
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4vr 2 7 (^(A + m 2 ) 1 a j )}, 



where we have introduced the new dual surface variable *<jj = *a + q*m\j\ 
which is spanned (/-times on the monopole current j: 5*(Tj = q*j. The 
flux of the unit charged magnetic monopole can be taken out by q strings 
carrying the unit flux. The mass of the gauge boson 9 is m = qy/^y/p. 

The BKT-representation (9) of the AHM partition function (6) can be 
also transformed to the dual representation using simple Gaussian inte- 
grations. We use two dual compact fields *B (vector field) and (scalar 
field) in order to represent the closeness properties of the currents Vj and 
*j, respectively. We also introduce two dual non-compact fields, *F (vec- 
tor field) and *G (rank-2 tensor field) in order to get a linear dependence, 
correspondingly, on the currents *<tj and *j under the exponential function: 

/oo roo rn rn 

V*F / V*G / V*B / V*^ 
-OO J — OO J — 7T J — 7T 

E ]T exp{-*/3(*G,(A + m 2 )*G) - * 7 (*F, (A + m 2 )*F) 

+i(*F,*a j )+i(*G,*j)+i(*B,6*a j - q*j) - i(*Z,5*j)} , (10) 

where 

* P = T&^ , * 7= i6^' 

Note that in this representation the integer variables Vj and *j are no 
more restricted by the closeness relations. Therefore we can use the Poisson 
summation formula with respect to these variables and integrate out the 
fields *F and *G. Finally, we obtain the dual field representation of the 
partition function (6): 

V*B / V*i E E 

exp{-*/? (d*5 + 2ir*u, (A + m 2 ) {d*B + 2tt*u)) (12) 
- *7 (d*^ + q*B + 2vr*«, (A + m 2 ) (d*^ + q*B + 2ir*v)) } , 

where *u and *v are the integer valued forms defined on the plaquettes and 
links of the dual lattice, respectively. Clearly, this is the dual Abelian Higgs 
model with the modified action. The gauge field *B is compact and the 
radial variable of the Higgs field is frozen. The model is in the London limit 
and the dynamical scalar variable is the phase of the Higgs field *£. 
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Thus in the presence of the dynamical matter the dual gauge field *B 
becomes compact 2 . The compactness of the dual gauge field implies that it 
is transforming under the gauge transformations of the following form: 

*B + d*a + 27T*k, (13) 

where the integer valued field k is chosen in such a way that *B £ (— n, it]. 

One can easily check that the operator (4) is not invariant under these 
gauge transformations: 

The invariance of the operator (4) under the gauge transformations (13) 
can be achieved if and only if the function *H X is an integer-valued form. 

Thus, if we take into account the Maxwell equation 5*H X = *5 X , we 
find that *H X should be a string attached to the monopole ("Mandelstam 
string"): *H X — > *j x , *j x € Z5, 5*j x = *6 X . The string must belong to the 
three-dimensional time-slice. However, one can show [7] that for a fixed 
string position the operator $ creates a state with an infinite energy. This 
difficulty may be bypassed [7] by summation over all possible positions of 
the Mandelstam strings with a measure ni*j): 

$ mo„,ncw = ^ £ l n \j,.) t ' ' r -' ! ' . (15) 

S*jx=*$x 

If Higgs field 4> is g-charged (q G 2Z), the summation in eq.(15) should be 
taken over q different strings with the unit flux. 

An example of a "reasonable" measure /J,(j x ) is [7]: 

^*j x ) = exp{-^\\*j x \\ 2 }. (16) 

This measure corresponds to the dual formulation of the 3D XY-model 
with the Villain action: 

S(x,r) = ^||dx-27rB + 27rr| | 2 . (17) 

Due to the compactness of the spin variables x t ne model (17) possesses 
vortex defects which enter the XY-partition function with measure (16). 
We thus defined "old" (4) and "new" (15) monopole creation operators. 

2 Another way to establish this fact is to realize that the pure compact gauge model is 
dual to the non-compact U(l) with matter fields (referred above as the (dual) Abelian 
Higgs model). Reading this relation backwards one can conclude that the presence of the 
matter field leads to the compactification of the dual gauge field *B. 
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3. Numerical results 

Below we present results of the numerical simulation of the new monopole 
creation operator. We investigate it in the simplest model which contains 
both the monopoles and the electrically charged fields: the Abelian Higgs 
model with compact gauge field and with the potential on the Higgs field 
corresponding to the London limit. The partition function for this model 
is given in eq.(6). The model has a nontrivial phase structure and we 
study both the phase where monopoles are condensed and the phase where 
monopoles are not condensed. 

First we substitute the monopole creation operator (15,16) into the dual 
representation of the compact AHM (12). Then we perform the transfor- 
mations back to the original representation: 



where we used the Wilson form of the action which is more suitable for 
the numerical simulations. The current j = *( 3 )*( 4 )j means that duality 
operation was first applied in the 3D time slice and then in the full 4D 
space. We have fixed the unitary gauge therefore the Higgs field was eaten 
up by the corresponding gauge transformation. 

The value of the monopole order parameter, ((f)), corresponds to the 
minimum of the (effective constraint) potential on the monopole field. This 
potential can be estimated as follows: 



We simulated the 4D Abelian Higgs model on the 4 4 , 6 4 , 8 4 lattices, for 
7 = 0.3. The larger the charge of the Higgs field, q, the easier the numerical 
calculation of I4fj( ( l > ) is. We performed our calculations for q = 7. For each 
configuration of 4D fields we simulated 3D model to get the ensemble of the 
Mandelstam strings with the weight ^(j x )- We generated 60 statistically 
independent 4D field configurations, and for each of these configurations 
we generated 40 configurations of 3D Mandelstam strings. We imposed the 
anti-periodic boundary conditions in the 31? space (the single monopole 
charge can not exist in the finite volume with periodic boundary condition). 

As we have stated above the weight function (16) corresponds to the 
3D XY-model with the Villain action. This model has the phase transition 
at k c (B = 0) « 0.32 [11]. Our numerical observation has shown that in 




(18) 




(19) 
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presence of the external field B, the critical coupling constant gets shifted: 
/ef « 0.42. 

Two configurations of the Mandelstam strings which correspond to con- 
densed (large k) and non-condensed (small k) phases of these strings are 
shown in Figures 1. 




Figure 1. Charged loops appearing in auxiliary theory for (a) non-condensed strings 
(k — 0.3) and (b) condensed strings (k — 0.5). 



We can expect that the operator (15) plays the role of the order pa- 
rameter in the n > k c phase, where the Mandelstam strings are condensed 
(k > k c ). In Figures 2 we present the effective potential (19) in the con- 
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Figure 2. The effective monopole potential (19) in (a) confinement and (b) deconfine- 
ment phases. 
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finement ((3 = 0.85) and deconfinement ((3 = 1.05) phases. The potential 
is shown for two values of the 3D coupling constants re > re c correspond- 
ing to high densities of the Mandelstam strings. In the confinement phase, 
Figure 2(a), the potential V(§) has a Higgs form signaling the monopole 
condensation. According to our numerical observations this statement does 
not depend on the lattice volume. In the deconfinement phase, Figure 2(b), 
the potential has minimum at <3? = which indicates the absence of the 
monopole condensate. 
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Figure 3. The effective monopole potential (19) in the 1ow-k region of the 3D model. 

For small values of the 3-D coupling constant re (in the phase where 
Mandelstam strings j x are not condensed), we found (Figure 3) that the 
potential V($) has the same behaviour for the both phases of 4D model. 
Thus the operator (15) serves as the order parameter for the deconfinement 
phase transition, if Mandelstam strings are condensed, i.e. re should be 
larger than k c {B). 

Summarizing, the new operator suggested in [7] can be used as a test of 
the monopole condensation in the theories with electrically charged matter 
fields. Our calculations indicate that the operator should be defined in 
the phase where the Mandelstam strings are condensed. The value of the 
monopole field which corresponds to the minimum of the effective potential, 
is zero in deconfinement phase and non zero in the confinement phase. 
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